Tailoring the flow of soft glasses by soft additives 
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We examine the vitrification and melting of asymmetric star polymers mixtures by combining 
rheological measurements with mode coupling theory. We identify two types of glassy states, a 
single glass, in which the small component is fluid in the glassy matrix of the big one and a double 
glass, in which both components are vitrified. Addition of small stax polymers leads to melting of 
both glasses and the melting curve has a non-monotonic dependence on the star-star size ratio. The 
phenomenon opens new ways for externally steering the rheological behavior of soft matter systems. 

PACS numbers: 82.70.-y, 64.70.Pf, 83.80.Hj 



The design of materials with well-defined rheological 
properties and the ability to alter these at wish, by tun- 
ing suitable control parameters of the system, are issues 
of central importance in today's soft matter research. 
Experimental findings, accumulating at a fast pace, call 
for the identification and profound understanding of the 
several underlying physical mechanisms that control the 
ability of soft materials to support stresses or flow under 
shear 0|. In several situations, e.g., in coating or pro- 
cessing applications, it is necessary to dramatically alter 
the viscoelastic properties of a material. One possibility 
to do so in a controlled way, is to exploit the phenomenon 
of dynamical arrest. Indeed, close to a liquid-glass tran- 
sition, a small variation of external parameters produces 
a spectacular change in the elastic properties of the ma- 
terial without significantly affecting its structure 

Star polymers have emerged as an ideal model sys- 
tem for exploring the flow behavior of soft matter and 
elucidating its molecular origin. They consist of / poly- 
mer chains covalently anchored onto a common center 
[3- Star polymer solutions are chemically simple, well- 
characterized and physically tunable in their softness, 
building thereby natural bridges between hard colloids 
and polymers. This property stems from the particular 
form of the effective, entropic interaction between star- 
polymer centers 0, ^ , which has a logarithmic /Yukawa 
form at small/large separations. Technologically, star 
polymers are important in several applications such as 
their use as viscosity modifiers in the oil industry JSj or 
their novel applications as drug-delivery agents |a, |3| • 

In this work, we show how one can gain control over 
the the rheological properties of soft, repulsive, glassy 
materials via the addition of a second, soft, repulsive 
component, which is chemically identical to the first but 



smaller in size. We consider a star-polymer glass per- 
turbed by smaller star-polymer additives. We observe 
melting of the big-star glass, induced by the small stars, 
and an unexpected non-monotonic dependence of the 
critical amount of additives, needed to melt the glass, 
on the big-small star size ratio. We also perform the cor- 
responding calculations of the melting line, within the 
binary mode coupling theory (MCT) framework ^ |^. 
The theoretical melting line is in agreement with the ex- 
perimental one, reproducing the qualitative trends of the 
latter. We demonstrate that the non-monotonic behavior 
arises from two different mechanisms by which the pres- 
ence of the additive significantly affects the rheological 
properties, depending on the size ratios between the com- 
ponents. The former rests on the fluidity of the smaller 
component in the single glass formed by the larger com- 
ponent. The latter stems from the mutual soft repulsions 
in the double glass, in which both components are vitri- 
fled. The melting of such a double glass represents a novel 
physical process, in which a glassy component is liquified 
through the addition of a second glassy component. 

Binary mixtures of big and small (1,4-polybutadiene)- 
stars were prepared in toluene. Three different types of 
big stars were employed, having an average functionality 
/i = 270 arms and arm molecular weight Ma ranging 
from 18 000 to 42 000g/mol 10]. Soft glasses were ob- 
tained at concentration ci/cl = 1.4, where cl denotes 
the big-star overlap concentration. Additional measure- 
ments were also performed with big stars of function- 
ality /i = 128 with arm molecular weight of about 
80000 g/mol, which are very similar to the 270-stars 
|lll |. We express the mixture composition by the val- 
ues piaf of the big and p2'Ji of the small stars, where 
Pi, i = 1,2 are the respective number densities and ai 
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FIG. 1: Experimental kinetic 'phase diagram' of binary star 
mixtures, indicating regions of liquid and glass for different 
concentrations of added small star p2 and size ratios q. Large 
stars of functionalities /i = 270 and 128 were used, with 
densities respectively piai = 0.345 and 0.412. Open symbols 
denote a liquid and closed ones a solid (glass) state. Only the 
data points closest to the melting lines are shown here. The 
lines through the data, passing above the full symbols and 
below the empty ones, are guides to the eye. 



are the corona diameters of the big stars appearing in 
the effective interactions employed in the theory, coin- 
ciding with the stars' hydrodynamic radii Rh,i, i = 1,2. 
Small stars with three different functionalities, /2 = 16, 
32, and 64, and molecular weights Ma between 1200 and 
80 000g/mol [Hlillll, were synthesized as well. Size 
ratios q = Rh,2/Rh.i = cr2/fi varied from 0.15 to 1. 
The mixture preparation protocol consisted of creating 
the big star glass at fixed number density picrj* = 0.345 
(ci/c* = 1.4) and then adding small stars with cer- 
tain g at a desired density p2, under conditions of very 
gentle and prolonged stirring. In this procedure, the 
glass was broken and the mixture was left to 'equili- 
brate' again. For the /i = 128-sample the fixed density 
was piCTi = 0.412. Dynamical rheological measurements 
(time, strain and frequency sweeps) were carried out in 
order to identify the state of the particular samples (solid 
or liquid behaviour). A strain-controlled rhcometer was 
utilized in the cone-and-plate geometry (25 mm diameter, 
0.04 rad cone angle) and dynamic frequency sweep tests 
were conducted in the range 100 — O.lrad/s at 20 °C in 
the linear viscoelastic regime. 

The experiments carried out provide evidence of U- 
shaped melting curves in the plane spanned by the size 
ratio q and the number density p2 of the small additives, 
as well as quantitative distinctions of two types of glasses 
depending on the value of q. The glassy samples are char- 
acterized by the typical virtually frequency-independent 
elastic moduh G' (350 Pa < G' < 800 Pa) and respec- 
tive weakly-frequency dependent viscous moduli G" that 
exhibit a broad minimum (^G^^^^^^O Pa <^ ^min ^ Pa) 
[l5| . For small concentration of small stars, the system 
exhibits, moduli G" > G" and G' ~ a;°, characteristic of 
solid, glass-like behavior 

[II 111 113 ■ As the density of 
the added small star increases, there is a dramatic change 



in the mixture's viscoelastic response. The glass melts, 
and a fluid results, whose moduli are weaker by orders 
of magnitude. A compilation of rheological results ob- 
tained from the available samples of different q and p2 
is presented in Fig. 1 for both studied functionalities of 
the large stars. In order to have a clear representation 
of all results in one plot, only the data closest to the 
glass-liquid boundaries are shown. This figure represents 
kinetic 'phase diagrams' of binary star mixtures, at fixed 
/i and pi and varying small-star characteristics in the 
(g, P2)-parameter space. The results suggest that a U- 
shaped melting line separates the fiuid (above) from the 
glass (below) states. This trend is general, since it is 
independent of the big stars used in the experiment. 

Our theoretical analysis is based on a coarse-grained 
star description employing effective interactions |5| in bi- 
nary star mixtures (181] with functionalities fi and corona 
diameters o'i.i = 1,2, combined with mode coupling the- 
ory (MCT) 0,0. In the pure star solution, MCT predicts 
a glass line at large / above an /-dependent density [31 , 
in agreement with experiment [2| • The equilibrium struc- 
ture factors Sij{k), k = 1,2, were calculated by solving 
the two-component Ornstein-Zernike equation |22l | within 
the Rogers- Young closure [2^ . The dynamics was calcu- 
lated using two-component MCT 0, . 

Upon addition of the smaller stars, the partial static 
structure factor S'ii(fc) of the large ones shows a loss 
of structure, a phenomenon caused by the depletion- 
induced softening of the repulsions between the big stars. 
At the same time, the small-stars partial structure fac- 
tor, £'22(^)7 gradually develops growing peaks. These 
structural changes result into melting of the glass at a q- 
and /2-dependent density Our theoretical results 

are summarized in Fig. 2, which can be interpreted as 
the theoretical 'kinetic phase diagram' of the system, in 
analogy with its experimental counterpart of Fig. 1. In 
agreement with experimental results, a characteristic Li- 
shape is found for the melting curves at all studied /2. 
Above q = 0.4 for /2 — 64, and q = 0.45 for /2 = 16 and 
32, no melting is found for any density of the small stars. 

The U-shape of the melting line points to the existence 
of two distinct microscopic melting mechanisms. For very 
small size of the additives, a standard depletion mecha- 
nism takes place. The osmotic pressure from the mobile, 
small stars, which are free to diffuse in the matrix formed 
by the big ones, leads to a reduction of the repulsive inter- 
actions between the latter. As a consequence, the cages 
that stabilize the glass are weakened, eventually breaking 
at a sufficiently high small-star concentration, at which 
the system melts. Since depletion is stronger at fixed p2 
as q grows, the melting curve has, in this regime, a nega- 
tive slope in the (g, /92)-plane, generating the left arm of 
the U-shaped line. However, as q further increases, the 
smaller stars become themselves slower and they begin 
to actively participate in the glass formation. At low ad- 
ditives concentration, the latter become trapped in the 




FIG. 2: Theoretical kinetic 'phase diagram' of binary star 
mixtures, calculated using MCT. The large-star functional- 
ity and concentration are fixed at the values /i = 270 and 
pi(Ti = 0.345. The diagram is shown for three different func- 
tionalities /2 of the small stars. Circles: /2 = 16; squares: 
/2 = 32; triangles: /2 = 64. The lines going through the 
calculated points are guides to the eye. The cartoons dis- 
play local arrangements in a single big-star glass with mobile 
small stars (left) and in a double glass in which there is mutual 
caging of both components (right). 



voids left out of the glassy matrix. Hence, two compet- 
ing mechanisms are at work: as before, the structure of 
the big stars is weakened by the addition of the small 
ones but, at the same time, the second component be- 
comes increasingly glassy. The onset of this mechanism 
brings about a reversal of the slope of the melting curve, 
since now the tendency of the small stars to soften the 
repulsion between the big ones is counter-driven by their 
own opposing tendency to jam. As p2 grows at fixed q, 
the small-star jamming cannot persist indefinitely, since 
the direct repulsions between the additives prevent them 
from occupying the same region of free space. As a result, 
melting of the glass takes place. We can distinguish be- 
tween two different glassy states under the U-curve, one 
in which only the big stars are jammed (low q) and one 
in which both components are arrested (large q) . We call 
the former single-glass and the latter double-glass state. 

The above interpretation is supported, first, by a com- 
parison with the results obtained if a one-component 
MCT-treatment is adopted, in which the smaller stars are 
assumed to form an ergodic fluid (dotted lines in Fig. 2). 
The two approaches yield very similar results at small size 
ratios, whereas the discrepancy becomes large at higher 
g's, signaling the tendency of the small component to ar- 
rest. The deviation between the two approaches can also 
be understood in terms of the scaling of the short-time 
mobilities with size ratio as q — > 1 the time-scale 

separation between the two species disappears. 

Further, we have calculated the elastic modulus G" fol- 



FIG. 3: Plot of the ratio G'/Gq, where Go is the modulus 
of the system without additives, against size ratio q close to 
melting. Theory refers to the case /2 = 32, while experiments 
refer to a compilation of several measured elastic moduli of 
star-star mixtures, obtained with different big and small stars 
at various concentrations of the additives. Inset: Dependence 
of the elastic modulus G' of the glass on p2 (/2 = 32). 



lowing Ref. 25]. In the inset of Fig. 3 we show the the- 
oretical results for G", demonstrating that, for q < 0.3, 
G' has a maximum at p2 = 0, whereas for the larger 
size ratio, q — 0.4, it has a minimum there. At small q, 
the additives lower G" through the softening of the big 
cages, whereas at high q they lead to stiffening of the 
glass through the fact that they are themselves driven 
to dynamical arrest. In Fig. 3, theoretical and exper- 
imental results for the normalized modulus G'/Gq are 
shown where Gq ^ 500 Pa stands for the average value 
of the elastic modulus of the big-star glass, without any 
additives. Experimental results display the same trends 
predicted theoretically. For small q, experimental values 
are much lower than the theoretical ones because MCT 
is not capable of taking into account the vast discrepancy 
in the mobilities of the two species in highly asymmetric 
mixtures. However, at larger g, in a full binary regime 
where the theory works best, the agreement between the- 
ory and experiments becomes almost quantitative. 

Additional evidence for the fluidity of the small stars 
at low g-values and their jammed nature at high ones 
is offered by the calculated nonergodicity factor /22(fc), 
shown in Fig. 4. Whereas the big-star nonergodicity fac- 
tor /ii(fc) is rather insensitive to q, and remains roughly 
the same in all glassy states, f22{k) shows a dramatic 
change. For q < 0.2, f22{k) is very small and its nonzero 
values are confined to a very narrow, small fc-domain. 
Such non-ergodicity factors are fully consistent with a 
mobile small component |23| . At high values 

/22(fc) has the typical range of a fully arrested system. 

The phenomenon of glass melting through polymer ad- 
dition has been observed in hard-colloid-polymer (CP) 
mixtures [2^ l3fil | . The physics and implications of this 
are, however, in our case very different. Whereas in the 
CP-case the solid melts due to the induction of short- 
range attractions on the colloids, which eventually drive 
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FIG. 4: The nonergodicity factors /ii(fc) (top) and /22(fe) 
(bottom) of the two components, for /2 = 32, p20"i = 0.02 at 
several q- values. For q = 0.3 the system is melted. 



the formation of a reentrant, attractive glass, in our case 
effective attractions are completely absent. This funda- 
mental difference is witnessed by the fact that in CP- 
mixtures melting takes place only for very asymmetric 
size ratios, q < 0.1, whereas in the present case it persists 
up to g ^ 1 . The melting of the solid in CP-mixtures re- 
quires the fluidity of the added polymer ji^l . To the con- 
trary, for star mixtures, the addition of a second com- 
ponent that becomes increasingly glassy is capable of 
bringing about a melting of the double glass. This is 
made possible by the crucial role played by the softness 
of the interactions between all species, which allows for 
the and rearrangement of soft cages, a mechanism absent 
in CP-mixtures. In star-linear mixtures, only the left arm 
of the U-curve is seen experimentally, the melting phe- 
nomenon ceasing altogether at g = 0.5 [IJ- There, the 
linear chains form a transient physical network, whereas 
stars are forced to maintain their spherical shape. 

Our findings are general, since they are based on 
the softness of the effective interactions involved. A 
large class of composite soft materials, such as charge- 
stabilized colloidal dispersions, microgels, dendritic poly- 
mers or self-organized micelles are characterized by soft 
effective interactions. Charged suspensions and dusty 
plasmas, for instance, can be described by Yukawa po- 
tentials whose decay length and strength can be tuned 
by controlling the charges and the screening microions, 
in a similar way that the strength and range of the ef- 
fective interactions between star polymers can be tuned 
by changing their functionality. Therefore, they should 
be amenable to manipulations of their rheology along the 
lines reported here. Functional versatility and viscosity 
control are desired material properties in medical, phar- 
maceutical and coating applications as well. Controlling 
these products near vitrification transitions, provides a 
means to achieve the desired properties variation. 
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